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Abstract. The effect of disorder is investigated in granular superconductive 
materials with strong and weak links. The transition is controlled by the interplay 
of the tunneling g and intragrain gintr conductances, which depend on the strength of 
the intergrain coupling. For g <C gintr the transition involves first the grain boundary, 
while for g ~ gintr the transition occurs into the whole grain. The different intergrain 
coupling is considered by modelling the superconducting material as a disordered 
network of Josephson junctions. Numerical simulations show that on increasing the 
disorder, the resistive transition occurs for lower temperatures and the curve broadens. 
These features are enhanced in disordered superconductors with strong links. The 
different behaviour is further checked by estimating the average network resistance 
for weak and strong links in the framework of the effective medium approximation 
theory. These results may be relevant to shed light on long standing puzzles as: (i) 
enhancement of the superconducting transition temperature of many metals in the 
granular states; (ii) suppression of superconductivity in homogeneously disordered films 
compared to standard granular systems close to the metal-insulator transition; (iii) 
enhanced degradation of superconductivity by doping and impurities in strongly linked 
materials, such as magnesium diboride, compared to weakly-linked superconductors, 
such as cupratcs. 
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1. Introduction 

The interplay of superconductivity and disorder has intrigued scientists for several 
decades pQ. Disorder is expected to enhance the electrical resistance, while 
superconductivity is associated with a zero-resistance state [2]. Bardeen, Cooper 
and Schrieffer explained the microscopic foundation of superconductivity in terms of 
pairing of electrons and the emergence of a many-body coherent macroscopic wave 
function [3]. Electron pairing defines a global order parameter A whose amplitude 
tends to zero by increasing temperature, current or magnetic field thus destroying 
the superconducting state. Anderson showed that weak disorder cannot lead to the 
destruction of the pair correlations. The transition temperature T c is insensitive to the 
elastic impurity scattering under the hypothesis that Coulomb interaction effects and 
mesoscopic fluctuations are negligible [H Ej. However, experiments performed on thin 
films have demonstrated a transition from the superconducting to insulating state with 
increasing disorder or magnetic field. In sufficiently disordered metals, these effects 
become important and the Anderson theorem is violated [61 [7J El [91 ITU] . 
Studies performed on homogeneously disordered conventional materials show, upon 
increasing disorder, the suppression of the superconducting critical temperature T c , 
the enhancement of the spatial fluctuations in A and the growth of the A/T c ratio 
[TT1 [T2| IT3] . More recently, impurity effects have been investigated in unconventional 
(i-wave superconductors, with the disorder causing pair breaking and suppression of 
T c [HI [151 HH [HI [THl EE]- The two-gap superconductivity is also affected by disorder. 
Experiments in neutron-irradiated MgB 2 show that the two-gap feature is evident in the 
temperature range above 21 K, while the single-gap superconductivity is well established 
as a bulk property, not associated with local disorder fluctuations, when T c is lowered 
to 11 and 8.7 K. The irradiation yields samples with extremely homogeneous defect 
structure so that the superconducting transition remains extremely sharp even in the 
heavily irradiated samples [201 IZE] ■ 

A still open issue in superconductivity is the enhancement of the critical transition 
temperature T c when some metals are in the granular forms rather than as a 
homogeneous bulk. It has been found that the enhancement is strongly dependent upon 
the intergrain coupling by varying pressure [221 G3], with many experiments confirming 
this phenomenon [2U [251 [261 [271 12H1 [291 E0] • Suppression of superconductivity in vicinity 
of the metal-insulator transition has been observed in homogeneous superconductors 
as amorphous Au x Sii_ x and Nb x Sil — x |31] . Chemical substitutions and impurities 
in MgB2 have resulted in superconductivity degradation and broadening of the R(T) 
curve pointing to an increasing effect of the disorder in such a strongly linked class of 
superconductors [32[33[Ml[3g[ia[33[I2 

Arrays of Josephson junctions with well controlled parameters are a very active 
field of research. As well as being of interest in their own right, they are also being used 
to model complex phenomena as a tool to investigate the effects of disorder in granular 

films [ISIlSEolEIlEalEaEl^ 
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This work is aimed at investigating the role of disorder in granular superconductors 
with different intergrain coupling, due to the presence of either strong or weak links. 
A parameter relevant to charge-carrier transport in such materials is the dimensionless 
tunneling conductance g = G/(e 2 /h), where G is the average tunneling conductance 
between adjacent grains and e 2 /h the quantum conductance. Films with g ^> 1 can be 
modeled as arrays of resistively shunted Josephson junctions, whose state is controlled 
only by the value of the normal resistance, rather than by the Josephson and Coulomb 
energies which are respectively defined as Ej = (jr/2)gA and E c = e 2 /C, with C the 
grain capacitance. The tunneling of normal electrons, which additionally takes place, 
results in the screening of the Coulomb energy, which reduces to the effective Coulomb 
energy E c = A/(2g). By comparing the Josephson energy to the effective Coulomb 
energy, one can notice that Ej is always larger than E c for g>l. This condition ensures 
the onset of the superconducting state at low temperature. Experiments show indeed 
that samples with the normal state conductance larger than the quantum conductance 
(i.e. with g 3> 1) always become superconducting at low temperature. 

A second parameter relevant to the understanding of the behavior of different 
granular materials is the intragrain conductance gi n t r - For standard granular systems, 
the condition g <C g intr holds. The intragrain region remains in the superconducting 
state, with the resistive transition occurring only at the grain boundaries. The condition 
g ~ gintr holds for tightly coupled grains, corresponding to homogenously disordered 
materials having comparable values of the bulk and grain boundary conductances 

[33J EH E21 EH E21 [63]. 

The different role played by the tunneling and intragrain conductances is 
determined by the strength of the coupling between the grains. In this paper, the 
conditions g <C g intr and g ~ g intr are considered in details. 

An array of Josephson junctions with different intergrain coupling and disorder 
degree is used to model the granular superconductor. The different contribution of g and 
gintr is accounted for by a proper circuital representation of the grain and its boundary 
within the network. The study is carried out by means of a numerical simulation whose 
main steps are summarized in Section [2j It is worthy of remark that the simulations 
reported in this work are carried out by the same numerical approach of Ref . [55] , where 
the different correlation shown by the current noise power spectra as a function of the 
intergrain coupling was investigated. The numerical results concerning the transition in 
weak- and strong-link networks as a function of the disorder are reported in Section |3] 
The transition temperature T c is lowered and the shape of the transition curve becomes 
smoother by increasing the disorder. Importantly, it is found that the disorder affects 
more dramatically networks with strong intergrain coupling. In Section HJ the results are 
quantitatively accounted for by estimating the resistive changes in weakly and strongly 
linked networks according to the effective medium approximation. 
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2. Numerical model 

As stated in the Introduction, the main purpose of this work is the investigation of 
the role of disorder in the resistive transition of granular superconductors with different 
intergrain coupling. The study will be carried out by adopting the numerical approach 
reported in Ref. [55], whose main steps are summarized here below. 

The resistive transition is simulated by solving a system of Kirchhoff equations for 
a network of nonlinear resistors biased by direct current, as shown in Figure E^a). Two 
types of networks are considered for describing the different intergrain coupling. The 
first type is the weak-link network for simulating materials, whose transition occurs in 
two subsequent stages. First, at low temperatures, the weak-links and, then at slightly 
higher temperatures, the whole grain undergoes the transition reaching the normal state. 
The weak-link network is used to model the first stage of the transition occurring at the 
grain boundary, while the grain interior still remains superconductive. The strong-link 
network is used for modelling the transition involving the whole grain. 

Grains are represented by a couple (triple) of nonlinear resistors for two-dimensional 
(three-dimensional) networks of Josephson junctions as shown respectively in Figures 



1 (b) and 1 (c) The nonlinear resistors give a basis of independent components of the 
current density able to reproduce the current flowing through the grain in arbitrary 
directions. The nonlinear resistors have current-voltage characteristics as shown in 
Figure |2] for underdamped (a), overdamped (b) and general (c) Josephson junctions. 
The Stewart-McCumber parameter (3 C = trc/ t j, where trc an d tj are respectively the 
capacitance and the Josephson time constants, identifies the three types: C ^> 1 (a), 
C <C 1 (b) and /3 C ~ 1 (c). The dependence of critical current I c ^j and magnetic field 
H C} ij on temperature can be written in the simplified form as: 



h,ij{T) =1 



co,ij 



H c ,ij(T) = H, 



CO,l] 



l-il 



Tr 



(la) 



(lb) 



where I C o,ij and H co ^j are respectively the low-temperature critical currents and magnetic 
fields and the exponent 7 ranges approximately from 1 to 2 depending on material 
properties. 

The current flowing through each nonlinear resistor defines the state (superconduc- 
tive, intermediate, normal) of the grain according to the current-voltage characteristics 
of the Josephson junction. As already stated, the disorder is introduced in the calcu- 
lations by random distribution of the critical current. The anisotropy is neglected and 
the same size is assumed for the grains. The reason for these simplifying assumptions is 
that these two features may additionally alter the network topology with a strong effect 
on the transition. In particular for small grain size, the values of the critical current 
might be correlated in neighboring grains. Therefore, the correlation length of disorder 
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should be taken into account by adopting a suitable spatial dependence of the critical 
current distribution. The critical current I Cji j(T) and the normal state resistance R ,ij 
are defined for each branch of the network. The intermediate state is characterized 
by the critical current I c>i j{T) and voltage drop between and V C)i j(T). The normal 
state, characterized by the resistance R ,ij, is reached when the current / crossing the 
Josephson junction exceeds I Cj ij. The disorder is introduced by taking the critical cur- 
rent I co ,ij & s a random variable distributed according to a Gaussian distribution with 



mean value I co and standard deviation ai c = y^2ij{Ico,ij — Ico) 2 /N '. Analogously, the 
disorder could be introduced by taking the critical field H coi j as a random variable, 
if the transition would be driven by an applied magnetic field H. The values of the 
resistances R^ between nodes % and j are taken as follows: 

Rij = if Vij ~ (superconducting state) (2a) 

Rij = Vij/I^ij if < < V C)i j (intermediate state) (2b) 
Rij = R Qy ij if > V Ct ij (normal state) , (2c) 

where V^ is the voltage drop between nodes % and j. The current- voltage characteristics 
is used to find the value of the voltage V^ and current I C) y by means of an iterative 
routine solving the Kirchhoff equations for the network. 

For weak-link networks, the resistance values R^ are calculated in a straightforward 
manner: the potential drops at the ends of each weak-link are compared to the potential 
values in the current- voltage characteristics according to Eqs. ff2a] I2b l I2c|) . Therefore, 
weak-links being respectively in the superconducting, normal or intermediate state can 
be distinguished. 

For strong-link networks, the resistance values R^ are calculated taking into account 
that the voltage drop across each grain is given by: 

1/2 



V, 



3 



(3) 



where corresponds to the voltage drop across each resistor Ru with j = 2 or j 



respectively for two- and three-dimensional arrays as shown in Figs. 1(b) and l(c^ 



Calculations are performed iteratively. First, a tentative set of potential values is 
chosen for all the nodes. Then, the resistance values Rij are calculated by using the 
Josephson junction current-voltage characteristics for any resistor between nodes i and 
j. Once the R^ are settled, the conductance matrix with entries Gy = 1/ Rij is defined 
and the new vector W± of the node potentials is calculated. The set of node potentials 
is introduced in the iterative routine and an updated vector W 2 is calculated. The it- 
eration is repeated until the quantity e n = \W n — W n -i\/\W n \ becomes smaller than a 
value e m i n chosen to exit from the loop. The simulations are performed by varying e min 
in the range 10~ 7 < e min < 1CT 11 to check that the value of e min does not appreciably 
change the results. The network resistance R is then obtained by W n (l)/I, where W n (l) 
is the potential drop at the electrodes. 
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3. Numerical results 

In this section, the results of the numerical simulations for different degrees of disorder 
are reported. It is shown that disorder affects at a different extent weak- and strong-link 
networks. 

At the beginning the network is entirely in the superconducting state (this condition 
is guaranteed by taking g ^> 1). Subsequently, the transition is made to occur through 
one of these processes: 

• the temperature is kept constant and the bias current (or the applied magnetic 
field) is varied. When the current exceeds the critical current I c ^j (or the 
magnetic field exceeds the critical field i7 c ,ij), the superconductive grain evolves 
to the intermediate and, then, to the normal state. 

• the bias current (or the magnetic field) is kept constant and the temperature is 
varied. A temperature increase causes a decrease of critical current 7 Cj jj according 
to Eq. ffTaj) (or of critical field H c ^ according to Eq. ( Tib]) ) and, ultimately, causes 
the transition of the grain to the intermediate and, then, to the normal state. 

As already stated, the disorder is modeled by assuming that the critical currents 
are a random variable distributed according to a Gaussian function with standard 
deviation aj c . The spread of the distribution of the critical currents determines the 
slope of the transition curve [61] . The standard deviation a Ic = corresponds to a fully 
ordered network, with all the Josephson junctions having the same critical current with 
the transition occurring simultaneously all through the network. When the disorder 
increases (cr/ c increases), the Josephson junctions have a wider spread of I C;i j and the 
network resistance changes more smoothly. 



Figures 3(a) and 3(b) show the resistive transition of the network for different values 
of o~i c for weak and strong link networks respectively. The temperature increases while 
the external current I is kept constant. As temperature increases, the critical current 
I Cj ij decreases according to Eq. (I lap . Links with I C)i j values smaller than 1^ undergo 
the transition to the normal state. If o~i c is small the resistive transition is steeper. 
In the limit of aj c = (no disorder in the network), the transition is vertical, since 
all the Josephson junctions become resistive for the same value of temperature. On 
the contrary, if o/ c is large the resistive transition broadens since the junctions become 
resistive at different temperatures. This effect occurs both in weak- and strong-link 
network, but is enhanced in strong-link networks. 



Figure 4(a) and 4(b) show the resistive transition when the bias current I increases 
at constant temperature, for different values of o~i c in weak- and strong-link networks re- 
spectively. When the bias current exceeds I c ,ij, the weak links become resistive. The 



transition curves of Fig. 4(a) and 4(b) exhibit a behavior similar to those of Fig. 3(a 



and 3(b) The disorder makes the resistive transition smoother, particularly in networks 



with strong-links. 
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4. Discussion 

In this section, the results of the simulations will be discussed. One can observe that 
the average network resistance R is determined by the elementary nonlinear resistances 
Rij between nodes % and j. The values of Rij depend on the external drive (current, 
magnetic field, temperature) and on the intrinsic properties of the junctions. The change 
of the resistance ARij can be expressed in terms of the external drive variation as: 

AR^^AI + ^AH+^AT . (4) 
The three terms on the right hand side of Eq. (jlj) can be written respectively as: 

^w AI - ~ t 1 ^ • < 5a > 
d Sr AH — H A ^> (5b) 

dRij „ m fdRijdlc dRa dH c \ A _ . . 

-w AT = [i>tw + diw) AT - < 5c > 

Equations ( l5al) and ( I5bl) mean that the increase (decrease) of bias current or magnetic 
field acts as a decrease (increase) of critical current I c or magnetic field H c . Eq. (l5cl) 
means that the temperature affects R^ mostly through a decrease of the critical current 
and magnetic field. By using equations (l5al l5bl l5cl) . with the derivatives dI c /dT and 
dH c /dT in Eq. (l5c|) calculated by using Eqs. fllalllbp . Eq. (jlj) can be rewritten as: 

AR « - -ff ( A '= + 4 AT ) - m ( AH < + ll k AT ) ■ (6) 

Equation ([6]) relates AR^ to the variation of critical current AI C or critical magnetic 
field AH C . One can note that ARij decreases when AI C or AH C increase due to the 
increased disorder in the array. Hence, since the network resistance R is proportional 
to terms varying as AR^, the slope of the resistive transition is smoother when AI C 
(AH C ) increases for a given temperature increase AT, regardless of the coupling strength 
between grains. 

However, Eq. cannot explain why the resistive transition becomes smoother 
with strong-link than with weak-links as one can notice in Fig. [3] and HI Therefore, in 
the following, the origin of the different behaviour exhibited by network with different 
intergrain coupling and same parameters of the elementary Josephson junctions, will be 
explained by including the effect of the different network topology. 

By effect of the temperature increase, layers of weak-links or grains either in the 
resistive or in the intermediate state, crossing the whole film are formed as shown is 



Figs. 5(a) and 5(b) The formation of a layer corresponds to an elementary step in the 
network resistance. This means that, in the limit of a large number of layers, which is 
a reasonable condition for real granular materials, the local slope AR of the transition 
curve can be approximated by the resistance of each layer Ri. In the remainder of this 
section, the resistance Ri will be estimate. Let N s j label the number of weak- or strong- 
links in the superconductive state before the transition of the layer. Let N a i label the 
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number of weak- or strong-links in the normal state and N m> i = N s j — N j the number of 
weak- or strong-links in the intermediate state at a given stage of the transition of each 
layer. The resistance Ri can be estimated as the parallel of the normal state resistors 
R ,ij and the intermediate state resistors R m ,ij a s ; 

R ,ijRm,ij / y \ 

1 ~ ~n~r Tn — r — ' ^ > 

The layer resistance Ri depends on the ratio of the normal N j and mixed state N m j 
resistances. For the strong links, the voltage drop between two neighboring grains is 
calculated according to Eq. (J3j) and thus is larger than Vy (voltage drop across each 
weak-link). Therefore, since the condition given by Eq. ( I2bl) is reached earlier, the 
denominator of Eq. (J7J) is larger in layers characterized by strong-links rather than 
weak-links for the same degree of disorder and bias current. This argument agrees with 
the fact that the resistive transition in strong-link networks occurs at temperatures lower 
than in weak-link networks. 



Fig. 6(a) shows the transition curves in weak- and strong-link networks with the 
same parameters. The slope is smaller for strong-link than for weak-link networks, 
consistently with the fact that the denominator of Eq. (J7J) is larger and thus AR w Ri 
is smaller. Furthermore, one can notice by comparing Figs. 6(b), 6(c) that the steps 
are higher for strong-links. This behavior has been confirmed by several runs of the 
transition simulations. Figs. [5] and [9] show nine samples of the resistive transition 
for weak and strong links respectively. One can clearly notice the different shape of 
the elementary steps. By implementing an automatic detection process of the steps 
endpoints, the elementary derivatives can be estimated. Fig. [TO] show the histograms of 
about 400 step slopes for weak-link (a) and strong-link (b) networks. This statistical 
analysis can be used for estimating an average value of the step slopes. The average 
ratio between derivatives for weak and strong links ranges between 1.3 and 2. A similar 
behavior is exhibited by the transition caused by current increase as shown in Fig. [7J 
To explain this issue, the elementary resistance between two neighboring sites i and j 
will be now estimated by using [65] [66], ETJ |68] : 

% = i?„exp(|^ + |), (8) 

where R Q = Tfcg/(e 2 7^),7^ is a rate constant related to the electron-phonon interaction 
lij ~ I); r ij 1S ^e distance between two sites, r Q is the scale over which the 
wave function decays outside the grain, is the zero field activation energy given by 
€ij = Ay(T) + E Ci ij, with E Cj ij = (3e 2 rij/(7re ed 2 ) the Coulomb energy and d the mean 
size of grain. Therefore, Eq. (jHJ) can be written as: 



A«(T) 



r, 



with 1/r* = [l/r + f3e 2 / \2ne ed 2 k B T)}. In Eq. (Q, the resistance R^ explicitly depends 
on the quantity r^, which is the effective distance seen by an electron flowing from 
grain % to j. The effective distance is different for electrons flowing either in weak- 
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or strong-link networks. Such a difference can be estimated by taking into account that 
at constant current the voltage drop V^ is proportional to r^. The voltage drop for the 
strong-link case is given by Eq. (J3j). A reduction of a factor Vij/[%2 V?] 1 / 2 of the distance 
Tij in comparison to the weak-link case should be correspondingly taken into account. In 
the simplest case of isotropic spherical grains, Vij is the same in any direction, thus the 
reduction factor is 1/ a/2 or 1/ a/3 respectively for two- and three-dimensional networks. 

By using the effective medium approximation [69], the average conductance G ema 
of the network can be calculated as follows: 

where z is the number of bonds at each node of the network, Gij = l/Rij and 
f(Gij) is the probability distribution function of the elementary conductance values 
G^. If the values G^ are continuously distributed according to the uniform function 
f(Gij) oc I /Gij, the average conductance is given by: 




The average conductance G ema varies as G2 times a factor depending to the ratio G\/Gi- 
The ratio G\/G<i is independent of the intergrain coupling contrarily to G<i- Therefore 
one can observe that the average conductance G ema increases as G2 increases with the 
coupling strength. According to presented model of the intergrain coupling, the value 
of the conductance G2 in case of strong and weak-link networks differs of the factor 
Kj/E ^ij] 1//2 - The average resistance R ema = 1/G ema is plotted in Fig. [TTJ It is worth 
noting that the resistance R em a for the strong-link case is always smaller then for the 
weak- link expected from the simulations. The presented discussion could be 

useful to explain existing experimental observations in granular materials that is very 
hard to understand with conventional mechanisms . (2U [25j [26j E3 EEJ EHJ 130] . 



5. Conclusions 



The effect of disorder has been studied in superconductors with different strength 
of the intergrain coupling. The superconductor has been modeled as an array of 
Josephson junctions, numerically solved by using Kirchhoff equations. The analysis 
shows that, on varying the external drive (temperature, current, magnetic field), the 
resistive transition occurs for lower T c and the R(T) curve broadens by increasing the 
disorder through a stepwise process. Importantly, it is found that the effect of disorder is 
more dramatic when the network simulates strongly rather than weakly coupled granular 
superconductors. The approach used and the results obtained in this work might add 
useful clues on the issue of the wide variability of critical temperature transition observed 
in real granular materials. It has been indeed observed an increase of critical temperature 
in compacted metallic powder compared to bulk samples of the same material. A 
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strong anticorrelation between the critical temperature T c enhancement and the value 
of metallic conductivity has been observed indicating that a major role is played by the 
electron-electron interaction which acts by suppression of the superconductivity [22j E3] . 
Chemical substitutions for Mg or B have been attempted to vary the superconducting 
transition temperature of MgB2- Most of the substitutions have produced a depression 
of T c and broadening of the R(T) curve contrarily to what observed in cuprates in 
which replacement of La by Y raises T c from 35K to 93K and sharpens the transition 
curve. It has been suggested that the two-band nature of MgB2 can result in an unusual 
behavior of its resistivity and T c as the material changes from the clean to dirty limits 
[321 1331 l34"l 1371 138] . The suppression/enhancement of T c is related to the competing 
effect of electron-electron and electron-phonon interaction which in their turn depend 
on size and radii of the compound and constituents. Intergrain and intragrain effects 
of disorder have been observed. Formation of magnesium or boron oxides result in 
poorly connected grains with an increase of intergrain resistivity and decrease of critical 
current density [351 [36J. At the same time, these oxides might migrate within the grains 
themselves increasing intragrain resistivity and flux pinning. Other impurities such 
as silicon, carbon, copper greatly affect critical current, temperature and resistivity 
[391 SB IH1 H21 H31 HH H5l 06]. Critical temperature degradation and broadening of the 
R(T)/R curve has been also observed in MgB 2 film by exposure to water [37]. The 
general feature of these experiments is that degradation of superconductivity seem to 
be related to the enhanced role of electron-electron interaction and impurity scattering 
in homogeneous metallic-like superconductors compared to the standard granular ones, 
i.e. that class of material whose intergranular conductance g is much smaller than the 
intragranular conductance gintv The dominant effect of the electron-electron interaction 
is taken into account in the present model by introducing a suitable circuital coupling 
among grains. 
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Figure 1. Scheme of a two-dimensional network representing the granular 
superconductor (a). Grains correspond to the nodes of the arrays. Each resistor 
of the network behaves as a Josephson junction with characteristics schematically 
shown in Fig. [2] For weak-link networks (YBCO-like materials), the grains remain 
in the superconducting state during the first stage of the transition. For strong-link 
networks (MgB2-like materials), the transition is a single-stage process involving the 
grains. The nonlinear resistors, characterizing each grain, are respectively shown for 
two dimensional (b) and three dimensional (c) case. 




Figure 2. Current-voltage characteristics for undcrdamped (a), ovcrdamped (b), 
general (c) Josephson junctions. For the general case (c), I m in,ij depends on the 
Stewart-McCumber parameter j3 c and ranges from and 0, where j3 c = trc/tj, 
where trc and tj are the capacitance and Josephson time constant respectively. 
p c » 1 (a), p c « 1 (b) and (3 C ~ 1 (c). 
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Figure 3. Resistive transition of two-dimensional networks with different degree of 
disorder at varying temperature for weak-links (a) and strong-links (b). The bias 
current / is kept constant. The degree of disorder is varied by changing the value of 
the standard deviation of the critical currents crj c from to 1 with step 0.1. 
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Figure 4. Resistive transition of a two-dimensional network with different degree 
of disorder at varying bias current for weak-links (a) and strong- links (b). The 
temperature T is kept constant. The degree of disorder is varied by changing the 
value of the standard deviation of the critical currents oi c from to 1 with step 0.1. 
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Figure 5. Scheme of two-dimensional networks when the first resistive layer is formed, 
respectively for weak- links (a) and strong- links (b). 
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Figure 6. Resistive transition of two-dimensional network with weak (blue) and 
strong (red) links as temperature increases. The standard deviation of the Gaussian 
distribution of the critical current aj a is equal to 0.2 for both curves. Zoom of the first 
step of the resistive transition in weak link (b) and in strong link (c) network. 



Resistive transition in disordered superconductors with varying intergrain coupling 19 




o 

2 




26.2 



26.4 



\l\co 

(b) 



a: 
c2 



0.025 
0.020 
0.015 
0.010 
0.005 
0.000 




AR=0.0032 



24.90 



25.05 



(c) 

Figure 7. Resistive transition of two-dimensional network with weak (blue) and strong 
(red) links as the bias current increases. The standard deviation of the Gaussian 
distribution of the critical current ui a is equal to 0.2 for both cases. Zoom of the first 
step of the resistive transition in weak link (b) and strong link (c) networks. 
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Figure 8. Resistive transitions of two-dimensional network with weak-links. 
Elementary resistance steps can be clearly observed. These nine curves are typical 
samples used for obtaining the data plotted in the histogram shown in Fig. 10(a) 
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Figure 9. Same as in Fig. [5] but for strongly linked grains. The histogram is shown 
in Fig. [lO(b)j 
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Figure 10. Histograms of the slopes of the elementary steps for weak-link (a) and 
strong-link (b) networks. The elementary steps have been obtained by means of 
transition curves similar to those shown in Figs. H] and El 
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Figure 11. Average network resistance R ema calculated according to the effective 
medium approach for z = 3 (a), z = 4 (b), z = 5 (c), z = 6 (d). Red curves refer to 
strongly coupled grains. Blue curves refer to weakly coupled grains. One can observe 
that the average resistance is smaller for strongly coupled networks for all the z values. 



